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Abstract 

We re-analyze data on deep inelastic polarized lepton-nucleon scattering, with particular 
attention to testing the Bjorken sum rule and estimating the quark contributions to the 
nucleon spin. Since only structure function data at fixed can be used to test sum 
rules, we use E142 asymmetry measurements and unpolarized structure function data to 
extract 5" at fixed = 2 GeV^. When higher-twist effects, which are important at low 
Q^, are included, both the E142 and SMC data are compatible with the Bjorken sum rule 
within one standard deviation. Assuming validity of the Bjorken sum rule, we estimate 
the quark contributions to the nucleon spin, finding that their total net contribution is 
small, with the strange quark contribution non-zero and negative. The quark spin content 
of the nucleon spin is in agreement with Skyrme model. 
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1 Introduction 



It was the Bjorken sum rule for polarized lepton-nucleon scattering, derived 

initially from quark current algebra at short distances, that led to the first pro- 
posal that deep inelastic structure functions should scale. However, although ap- 
proximate scaling has been observed in deep inelastic scattering and quark model 
relations for structure functions hold in general, with deviations which are under- 
stood quantitatively in terms of perturbative QCD and asymptotic freedom at large 
momentum transfers, the Bjorken sum rule was never actually tested during the 25 
years following its discovery. Data on polarized lepton-proton scattering have been 
available for about 15 years 0-05 but data on polarized lepton-neutron scattering 
have become available only very recently [§[,[§]. The stimulus for performing these 
experiments, finally, was provided by the results of the EMC experiment 01,0 on 
polarized muon-proton scattering, which differed from predictions based on naive 



constituent quark model ideas [jl^], but were consistent with perturbative QCD. 
Polarized lepton-neutron scattering would also have to differ from naive constituent 
quark model ideas if the Bjorken sum rule were to be valid. If it were not valid, the 
whole edifice of perturbative QCD and asymptotic freedom would collapse [jlT 



One of the two results on polarized lepton-neutron scattering announced recently, 
the one by the SMC |Q, does indeed find a deviation from naive constituent quark 
model ideas, and, when combined with the EMC proton data, verifies the Bjorken 
sum rule within admittedly large errors. However, the other polarized neutron 
experiment by the E142 collaboration |Q, which has much smaller statistical errors, 
is claimed, when combined with the EMC data, to disagree with the Bjorken sum rule 
by two standard deviations. This result therefore presents a puzzle whose resolution 
apparently needs affirmative answers to one or more of the following questions: are 
the EMC data wrong? are the E142 data wrong? is the Bjorken sum rule wrong? 



We take the point of view in this paper that one must be prudent before jumping 
to any such dramatic conclusion, especially before concluding that the Bjorken sum 
rule is wrong. We recall that in every polarization experiment the data have been 
taken in different ranges of Q^, which varied from bin to bin in the Bjorken scaling 
variable x, with average values that are known in the case of the SMC and 
EMC experiments, but have not been published by the E142 collaboration. The 
EMC and the SMC have used the fact, observed also in other experiments, that 
the asymmetry Ai{x) has no detectable dependence on in the range studied, 
to estimate the polarized structure functions at a fixed value of Q^. Whereas the 
EMC used their own measurements of the unpolarized structure function F2{x, Q^) 
and a perturbative QCD model for the ratio R{x, Q^) of longitudinal to transverse 
virtual photon cross-sections to calculate gi{x,Q'^) for Q'^ fixed at 10.7 GeV^, the 
SMC used the more recent NMC parametrization [|1^] of the unpolarized structure 
function F2{x, Q"^) and a SLAC parametrization [0 of the ratio R{x, Q"^) to calculate 
gi{x, Q"^) for fixed at 4.6 GeV^. This procedure is a priori better suited for testing 



1 



sum rule predictions, which are formulated for structure functions integrated over 
all X values at fixed Q^, than would be an integral over x at variable values of Q^. 



In this paper, we first check whether there is indeed any evidence against the 
validity of the Bjorken sum rule. To do this, we start by rescaling the EMC polarized 
proton data to at 4.6 GeV^ for comparison with the SMC H], and at 2 GeV^, 
for comparison with E142 P] using the latest NMC parametrization of the 
unpolarized structure function and the latest SLAC parametrization of R{x, Q'^) 
[|1|]. For the latter comparison, we use the measurements of reported by E142 
and the same NMC and SLAC parametrizations to extract gi{x, Q"^) for fixed 
at 2 GeV^. Next, we discuss the uncertainties in the extrapolations of the E142 
data to X = and 1. Since the appropriate power behaviour at small x is not 
tightly constrained by other data or by theory, this introduces some extra uncertainty 
beyond that discussed by E142. A priori, the error in the SMC extrapolation to x = 
is 5 times smaller, since they have measurements down to x = 0.006, as opposed 
to the 0.03 of E142. Contrary to what is stated in ref. §], the behaviour of the 
polarization asymmetry at large x is not determined by perturbative QCD without 
extra assumptions about the non-perturbative structure of the nucleon, which also 
tends to increase the uncertainty beyond that assumed by E142. Next we show that, 
when proper attention is paid to the dependence of the structure functions at 
large x, the SMC data are largely compatible with a recently derived bound M . 
We then extract the integrals of gf - Qi at fixed = 4.6 GeV by combining the 
rescaled EMC and SMC data, and at 2 GeV^ by combining the rescaled EMC and 
E142 data, including realistic allowances for the errors. Particularly in the low-Q^ 
case of the comparison of E142 and the EMC, higher-twist and mass corrections 
must be taken into account. According to the indicative estimate in ref. |T5| these 
corrections are more than 10% at = 2 GeV^ for the Bjorken sum rule. Q The 
estimates of ref. should therefore be included in the theoretical prediction. After 



including these subasymptotic effects, both the E142 as well as the SMC data are 
in good agreement with the Bjorken prediction. Encouraged by this agreement, 
we assume that the Bjorken sum rule is indeed correct, and go on to extract the 
contributions of the various quark flavours to the proton spin using all the available 
polarized structure function data. We find 

A?i + Arf + As = 0.22 ± 0.10 (1) 

which is far from the naive quark model. It is close to the prediction of the Skyrme 
model for light quarks Jl^. We conclude that a quantitative understanding of axial 



current matrix elements is now emerging. 



^See however ref. |16 for a recent discussion of higher-twist effects in this context. 
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2 Data at Fixed 



We first discuss the rescaling of tlie available polarized structure function data 
to a fixed value of Q^, which is a crucial step for testing QCD sum rules, that are 
formulated at fixed . Experiments measure directly the polarization asymmetry 
Ax{x^ Q"^), in terms of which the polarized structure function gi{x, Q"^) is given by 



Ai(x,g^)F2(x,Q^) 
2x[l + i?(x,Q2)] 



(2) 



where F2(x,Q^) is the conventional unpolarized structure function, and Q^) 
is the ratio of virtual longitudinal to transverse virtual photon cross sections. As 
already mentioned, each experiment takes data over a range of that varies with 
X, and differs from experiment to experiment: (Q^) = 2, 4.6, 10.7 GeV^ for E142, 
the SMC and the EMC respectively. Each experiment reports that the asymmetry 
that it measures exhibits no detectable variation with in any x bin. It is therefore 
reasonable to use this Q^-independent value, together with the unpolarized -F2(x, Q^) 
and R{x, Q^) measured in other experiments, to estimate gi{x, Q"^) at some fixed Q^, 
taken ideally to be the average (Q^) for each experiment. However, to combine the 
experiments that have different {Q"^) to test the Bjorken sum rule, it is necessary 
to estimate gi{x,Q'^) for at least one non-ideal value of Q^: the variation of 
gi^x.Q"^) is substantial, even if that of Q^) is negligible. 



Both the EMC and the SMC |§| have presented their data at fixed as 
advocated above. However, the EMC used a parametrization of the unpolarized 
data that has now been superseded by more recent results from the NMC and 
SLAC The latter have been used by the SMC to estimate their gf{x,Q'^), 

and we use them here to re-estimate gi{x,Q'^). Figure 1 shows the rescaled EMC 
gi{x, Q"^) at fixed = 2, 4.6 and 10.7 GeV^, using a smooth parametrization of 
the Ai EMC data. The two former values of are needed for combining with the 
E142 and SMC data respectively: the Q^-dependence of gi{x, Q^) is seen clearly, 
and the corresponding integrals Ti{Q'^) = Jq gi{x, Q^) are shown in Table 1. 



Table 1 



Q'(GeV^) 




2.0 


0.124 ± 0.013 ± 0.019 


4.6 


0.125 ± 0.013 ± 0.019 


10.7 


0.128 ± 0.013 ± 0.019 



The integrals r{{Q'^) = gl{x,Q'^) obtained from the EMC data on Al{x,Q'^), at 
Q2 = 2 GeV^, combined with the NMC data for F^{x, Q^) and SLAC parametrization 
of R{x^ Q^) = (Tl/'^t- The first is the statistical and the second is the systematic error. 
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We assume here that the systematic error on is as quoted by the EMC [0], even 
though the systematic error on the NMC Ff' is smaller than the systematic error on 
the EMC . If there were no error at all on the systematic error in would 
be reduced from ±0.019 down to ±0.018. The difference from the value quoted by 
the EMC for = 10.7 QeV^ is well within their quoted errors: we will comment 
later on the Q^-dependence of Ti{Q'^). 

We have also applied the same procedure to the E142 data , using their values of 
A^lx) and the NMC and SLAC parametrizations of F2{x, Q"^) and R{x, Q"^) respec- 
tively to estimate gi{x, Q'^) at Q'^ = 2 GeV^, their average value. This procedure 
yields 

/ dxg'lix, g2=2 GeV^) = -0.022 ± 0.006 (stat.) ± 0.006 (syst.) (3) 

J 0.03 

to be compared with the estimate of —0.019 ± 0.006 (stat.) ± 0.006 (syst.) given 
by E142. 



Before estimating F"(Q^ = GeV^) on the basis of equation (^), we first comment 
on the extrapolations of their data beyond x = 0.6 and below x = 0.03. Models for 
the polarization asymmetry that combine a non-perturbative Ansatz for the neutron 
wave function with perturbative QCD at large can be used to estimate the limit- 
ing value of A^{x) as x — > 1, but perturbative QCD alone does not predict a limiting 
value. Therefore, we prefer to estimate the error due to the high- a; extrapolation by 
allowing |y4^(x)| < 1, rather than by specifying a limiting value. Thus we estimate 
a contribution 0.000 ± 0.003 from the x > 0.6 region. The low-x extrapolation of 
the E142 data is a priori more uncertain than that of the SMC, because the latter 
measure down to lower x = 0.006. E142 have chosen to extrapolate using a power 



form gi{x) = Ax°', where a was fixed at the value 0.18 predicted in ref. [18| and 



adopted in ref. |T^. The model used in ref. [|IU| is however incompatible with the 



new data, and hence cannot be used to support the value of a adopted by E142. 
Other data do not fix the power a: for example, a fit to the EMC data at low x 
gave a = 0.07to;42 , whereas we consider ||21| a plausible theoretical range to be 



< a < 0.5, as assumed by the SMC. Using this latter range of a to estimate the 
contribution of the range x < 0.03 as —0.006 ± 0.006 0, we estimate the full integral 
to be 

F^(g2=2 GeV^) = -0.028 ± 0.006 (stat.) ± 0.009 (syst.) (4) 
to be compared with the E142 estimate of —0.022 ± 0.011. 



The value of TKQ^) at = 4.6 GeV^ can be extracted from the SMC data g] 
together with the appropriate value of F]'(Q^) in Table 1. This results in 

F^(Q2=4.6 GeV^) = -0.076 ± 0.046 (stat.) ± 0.037 (syst.) (5) 



^We note that one would obtain a lower value if one combined E142 and SMC data, which are 
lower at small x. 
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to be compared with the SMC estimate of -0.08 ± 0.04 ± 0.04. 



There has been some question [T^ whether the SMC data at large x are compatible 
with general positivity bounds derived within the framework of the quark-parton 
model. An interesting bound is derived in ref. [jl4| by eliminating the u quark- 
parton contributions to gi{x, Q'^) and Q'^)\ 



< ^d{x,Q'')+^s{x,Q^), 

(6) 

It should be noted that this bound could be violated by higher-twist effects, that are 
expected to grow at smaller and larger x. An additional possible source of error 
is the assumption, made in ref. that gi{x, Q^) -|- Q"^) = gf{x, Q"^) locally 



in X. The superficially similar relation among the first moments used by the SMC 
[|] r?(g2) + r^(Q2) ~ rf(Q2)/(i _ LScj^), where ud = 0.058 is the probability of 
the deuteron to be in a Z)-state, is well justified, but the local relation is violated 
by smearing of the deuteron at large x. 



We re-express the bound in terms of the directly-measured quantities v4f (x) 
and and the Q^-dependent quantities R{x,Q'^) and ^(x, Q^) = 

F^ix,Q')/Fi{x,Q'y. 



AAUx,Q'nx,Q')-Al{x,Q' 



< 



i + i?(x,g2) 4e(x,g2)_i 



(7) 



We have evaluated both sides of this version of the bound, using the same fixed 
values = 2, 4.6 GeV^ on each side to compare EMC data with the E142 and 
SMC data respectively. 

The E142 data satisfy the bound comfortably, whereas the SMC data are in 
marginal disagreement, as shown in figure 2. We do not take seriously the latter 
disagreement, in view of the large errors and the possible contributions of higher 
twist effects. Shifting the central values of the SMC data at x > 0.4 so as to be 
consistent with the bound would not in any case change significantly our estimate 
of r^(Q2) based on the SMC data. 



3 Testing the Bjorken Sum Rule 

We are now in a position to evaluate the difference between rf(Q^) and r"(Q^) 
at = 2 GeV^ using the EMC and E142 data, and at = 4.6 GeV^ using the 
EMC and SMC data: 

r^'-"(g2=2.0 GeV^) = 0.152 ± 0.014 ± 0.021 

(8) 

rP-"(Q2=4.6GeV2) = 0.201 ± 0.048 ± 0.042 
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These are to be compared with the Bjorken sum rule, which takes the following form 
when the leading-order perturbative QCD corrections are included 0: 



dx 



gf{x,Q')-g^{x,Q' 
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9a 



l-a,(Q2)/7r 



(9) 



To estimate the appropriate values of as{Q'^) at = 2, 4.6 GeV^, we take the latest 
estimate at ^^(m^) = 0.330 ± 0.046 |2^, which is consistent with determinations 
at higher Q^, and use the leading order corrections to calculate 0^(2, 4.6 GeV^) = 
0.371 ±0.058, 0.304 ±0.039, together with gA = 1.2573 ± 0.0028 p^, leading to the 
following values of the right-hand-side of equation @: 



rf-"(g2=2.0GeV2) = 0.185 ± 0.004 



rf-"(Q2=4.6GeV2) = 0.189 ± 0.003 



(10) 



However, these values cannot to be compared with the evaluations (|^) before in- 
cluding subasymptotic 0{1/Q^) effects due to mass corrections and higher-twist 
operators, which are particularly important at the low average of E142. 



The mass corrections ||T5[ are proportional to moments of gi'"'{x), which can be 
estimated using the data themselves: 



2\P 



dx x'^gi{x, Q^) 



0.0168, (Q2 = 2.0GeV^) 
0.0130, (g2 = 4.6GeV2) 



'11^ 



2\n\ 



dxx gi{x, Q 



< 7 X 10-^ 



Errors in {x^)p can be neglected, since the mass corrections include an additional 
small factor (4/9)(m^/(5^) where n^v is the nucleon mass. Estimating the higher- 
twist corrections requires calculating the reduced matrix elements {{U^^)) and 



{{U )) of the local operators of spin one and twist four [|23| , |[T5[] : 



^ ugG^^Yu + dgG^^Yd + — sgG^^Ys ■ 



;i2) 



[/^^^ = ugG^^ru-dgG^^rd., 



{N\U,\N) = s,{{U)), s,^N^,^,N. (13) 

where G^j, = e^uapG'^'^ Xa/2 and N, N are nucleon spinors. The matrix elements of 
the non-strange quark operators have been estimated using sum rule techniques in 
ISj]: the possible error quoted is estimated at 50% for rj*"" and 100% for Fj''''", and 



6 



it has been suggested that the matrix element of the strange quark operator should 
not be larger than this uncertainty. Ref. [|I5| found the following magnitudes for 
the higher-twist correction: 

(5r?""(Q2=2.0 GeV^) = -0.022 ± 0.011 
5rr"((5^=4.6 GeV^) = -0.010 ± 0.005 

(14) 

5rr"(<5^=2.0GeV2) = +0.017 ± 0.017 
5r?+"(Q2=4.6GeV2) = +0.007 ± 0.007 

The net contribution of the higher-twist and mass corrections to r^{Q'^) is quite 
small, becoming negligible at the of the EMC experiment, and has the effect 
of decreasing the integral slightly at smaller Q^. This is consistent with the trend 
found in Table 1 above on the hypothesis that is independent of Q^, though 
the inferred variation cannot be considered significant. The net contribution of the 
l/Q"^ corrections to T^IQ^) are much larger, positive, and particularly significant 
for the low of E142. Thus it is not surprising that their value of r"(Q^) happens 
to be larger than that of the SMC. The net correction to the Bjorken sum rule is 
important, exceeding 10% at Q"^ = 2 GeV^: 

rr"(Q2=2.0GeV^) = 0.163 ± 0.012 

(15) 

r?""(g2=4.6GeV2) = 0.180 ± 0.006 

Comparing these theoretical estimates with the experimental evaluations (^, we 
conclude that the EI42 and EMC data are consistent with the Bjorken sum rule 
within less than one standard deviation. This point is shown graphically in figure 3. 
One way of phrasing this consistency is to quote the effective values of qa extracted 
from the two combinations of experiments by subtracting the subasymptotic correc- 
tions (0) and removing the asiQ"^) correction appearing in (|^): 



EMC & SMC: gf = 1.39 ± 0.38 



EMC & E142: gf = 1.19 ± 0.17 



(16) 
(17) 



Combining all three experiments, we find 



gf = 1.22 ± 0.15 (18) 

The Bjorken sum rule is consistent with all the available data on polarized structure 
functions, and is now verified at the 12% level. 



7 



4 Evaluation of the Quark Contributions to the 
Nucleon Spin 

Having reassured ourselves that the available data on polarized nucleon structure 
functions are consistent with the Bjorken sum rule, we now extract the contributions 
of the u, d and s quark flavours to the total nucleon spin, and compare them with 
theoretical estimates. We treat independently the EMC data on polarized protons, 
the SMC data on polarized deuterons, and the E142 data on polarized ^He. In the 
case of the SMC result, we fit the data after applying the D-wave correction. In the 
case of the E142 result, we assume the same model of ^He as in ref. (see also 
ref. |2^ for a detailed discussion of ^He wavefunction in this context). Before fitting. 



we apply to each value of the integral of gi{x, Q"^) the appropriate higher-twist, mass 
and perturbative QCD corrections. We denote the resulting corrected first moments 
rf((5^) and r"((5^), which correspond to the net charge-weighted sums of quark 
helicities. Thus we have 

^ ^-AM + -Arf+-As) = f{ (g2=io.7GeV2) = +0.140 ± 0.023 



2 V9 9 9 
1 /I . 4 . , 1 . ^ ^„ 



-AM + -Arf + -As = (g'= 2.0 GeV') = -0.056 ± 0.015 (19) 

2 \9 9 9 / 

]- f ^Am + + ^ As) =fl + fl {Q^= 4.6 GeV^) = +0.048 ± 0.055 
2 \9 9 9 / 

We combine each of these inputs with 



An- Ad = = 1.2573 + 0.0028 (20) 



and take [p3] 



F = 0.46 + 0.01; D = 0.79 + 0.01; F/D = 0.58 + 0.02 (21) 

Each of eqs. (0) and (^) can be rewritten as linear constraint relating the values of 
AS and As. In fig. 4 we plot the allowed regions in the AE - As plane, correspond- 
ing to these constraints. The three data sets above yield the following compatible 
estimates of the total quark contribution to the proton spin: 

EMC (g2 = 10.7 GeV^): An + Arf + As = +0.16 + 0.21 

An + Arf + As = +0.06 + 0.25 (22) 
An + Ac/ + As = +0.29 + 0.14 

Ignoring any possible dependence in An + A(i + As, we combine these in quadra- 
ture to obtain 

AE = An + Arf + As = 0.22 + 0.10 (23) 
where the individual quark contributions are 

An = +0.80 + 0.04 

M = -0.46 + 0.04 (24) 
As = -0.13 + 0.04 



SMC (g2 = 4.6 GeV^) 
E142 (g2 = 2.0 GeV^) 
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The global estimate of As is approximately three standard deviations away from 
zero, providing good evidence that the naive assumption of ref. |]T0[, namely that 



strange quarks should not contribute to the nucleon spin, is false. This does not 
strike us as implausible, since we have much more sophisticated models of nucleon 
structure then when [jlOi was written almost 20 years ago. 



In particular, we note that the value (|23|) of SAg is quite close to zero, which is 
the value expected in the class of chiral soliton models pioneered by the computa- 
tion in the Skyrme model in the limit of massless quarks and a large number of 



colours. It has been suggested that corrections to this limiting value might amount 
to about 30% in the case of realistic quark masses and three colours which would 
be consistent with the experimental value (|23|). We are not aware of any other 
dynamical model of nucleon structure which predicts a value of SAg close to zero, 
though some models are able to accommodate it. For example, it has been sug- 
gested [^] that gluons inside the nucleon might be highly polarized, and that this 



polarization would communicate itself to the quarks via radiative corrections and 
the axial anomaly. Unfortunately, this interesting suggestion does not give insight 
why EAg should vanish in any limitQ. There have been some initial calculations of 
matrix elements "from first principles" using lattice techniques [^-[Q, which are 
consistent within errors with our determinations. 

It is also worthwhile noticing that As being nonzero is a part of an interesting 



pattern |^^: experiment indicates that certain strange-quark bilinear operators, 
such as S7^75S have relatively large matrix elements in the proton, while others are 
very small. The presence of a substantial non-valence component of ss pairs in the 
proton has some striking consequences. One of these is the evasion of the OZI rule 
in the couplings of ss mesons to baryons , leading to surprisingly large branching 



ratios for production in pp annihilation at rest [36 . 



5 Interpretation 

We conclude that experiment and theory are converging on consistent and quite 
quantitative decompositions of the spin of the nucleon. The chiral soliton models 
pioneered by the Skyrme model are perhaps unique in providing understanding of 
these results, which disagree with naive quark model ideas. The basic idea of these 
chiral models is that a nucleon contains a very large number of very light, relativistic 
quarks, that are best described by a topological lump in an effective bosonic field. 
The data seem to indicate that this is a better starting-point for understanding the 

•^See ref. ||2^ for discussion of \/Nc corrections and ref. |^ for a current assessment of the 
situation in Skyrme type models. 

"^For proposals in this direction, see ref. po|| . 
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spin of the nucleoli than the naive picture of three non-relativistic valence constituent 
quarks. However, the constituent quark model does provide good understanding of 
many other aspects of nucleon structure. We are therefore confronted with the 
challenge |^| of relating constituent quarks to the chiral picture that works so well 
for the polarized structure function data re-analyzed in this paper. There is already 
some indication [R3]-Pl[] that this program has good chances of success. 
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Figure Captions: 



Fig. 1. (a) The polarized structure function gi{x, Q"^) at = 2.0, 4.6 and 10.7 
GeV , obtained from the EMC data 0, smoothed and combined with the NMC 
[|l^ data for F2{x,Q'^) and SLAC parametrization |13| of R{x,Q^) = ai/aT- Con- 
tinuous curve: = 2 GeV^ dots: = 4.6 GeV^ dot-dash: = 10.7 GeV^ 
Compensation between the decrease at large x and the increase at small x results 
in the relatively small Q^-variation of TilQ"^) found in Table 1. (b) The polarized 
structure function g^{x,Q^) at = 2.0, 4.6 and 10.7 GeV^, obtained from the 
E142 data 0, smoothed and combined with the NMC data for F^{x,Q^) 
and FKx,^^), and SLAC parametrization [|T^ of R{x,Q^) = ai/aT- Continuous 
curve: = 2 GeV^ dots: = 4.6 GeV^ dot-dash: = 10.7 GeVl Unlike in 
the case of gi{x,Q'^), there is no compensation between the large and small x, and 
this results in the relatively large Q^-variation of r"(Q^). 



Fig. 2. The difference between the right-hand and left-hand sides of the bound (0). 
The actual errors are slightly larger than those indicated by the error bars, as the 
latter refer to the error in the left-hand side only, (a) E142 data , combined with 
EMC data §, rescaled to = 2 GeV^; (b) SMC data [|, combined with EMC 
data 0, rescaled to Q"^ = A.Q GeV^ 



Fig. 3. Experimental tests at = 2 GeV^: E142 and EMC, = 4.6 GeV^: SMC 
and EMC of the Bjorken sum rule, including perturbative QCD corrections (dot- 
dashed lines) and higher-twist corrections (solid lines). The asymptotic value qa/^ 
is denoted by a dotted line. 

Fig. 4. The allowed regions in the AS-As plane, corresponding to the linear con- 
straints (0) and (0). Continuous lines: AS - 3 As = 3F - D; dots: t^{Q^) 
constraint; dot-dash: f^(Q2); dashes: f{{Q^) +f'^{Q^). 
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